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ABSTRACT 

We have presented the non-equilibrium (time-dependent) cooling rate and ionization 
state calculations for a gas behind shock waves with v ~ 50 — 150 km s -1 (T s ~ 

0. 5 — 6 x 10 5 K). Such shock waves do not lead to the radiative precursor formation, 

1. e. the thermal evolution of a gas behind the shock waves are controlled by collisions 
only. We have found that the cooling rate in a gas behind the shock waves with 
v ~ 50 — 120 km s _1 (T s ~ 0.5 — 3 x 10 5 K) differs considerably from the cooling rate 
for a gas cooled from T = 10 8 K. It is well-known that a gas cooled from T = 10 s K 
is thermally unstable for isobaric and isochoric perturbations at T ^ 2 x 10 4 K. We 
have studied the thermal instability in a collisionally controlled gas for shock waves 
with v ~ 50 — 150 km s . We have found that the temperature range, where the 
postshock gas is thermally unstable, is significantly modified and depends on both 
gas metallicity and ionic composition of a gas before shock wave. For Z £ O.IZq the 
temperature range, where the thermal instability criterion for isochoric perturbations 
is not fulfilled, widens in comparison with that for a gas cooled from T = 10 s K, 
while that for isobaric perturbations remains almost without a change. For Z ~ Zq 
a gas behind shock waves with v ^ 65 km s _1 (T s £ 10 5 K) is thermally stable to 
isochoric perturbations during full its evolution. We have shown that the transition 
from isobaric to isochoric cooling for a gas with Z £ O.IZq behind shock waves with 
T s = 0.5 — 3 x 10 5 K proceeds at lower column density layer behind a shock wave than 
that for a gas cooled from T = 10 8 K. The ionic states in a gas with Z ~ 10~ 3 — 1 Zq 
behind shock waves with T s £ 4x 10 5 K demostrate a significant difference from these 
in a gas cooled from T = 10 8 K. Such difference is thought to be important for correct 
interpretation of observational data, but hardly help to dicriminate thermally stable 
gas. 



Key words: cosmology: theory - intergalactic medium - quasars: general - absorp- 
tion lines - physical data and processes: atomic processes 



1 INTRODUCTION 



T herma l instability (TI) comprehensively analyzed by 

iFieldl 1 19651 ) is frequently considered as a good candidate 
to ex plain the planet ary nebulae e volution ( Hunter fc Sofial 



Il97ll ), cooling fl ows dNulserJll986T). the formation of the in- 
terstellar clouds (|Burkert fc Linll2000l ). the temperature dis- 
tribution of the uns table interstellar gas in the Galactic disk 
jGazol et al. Il200ll). the formation of high velocity clouds in 
galactic corona " pinnev et al. 1 120091 ) and so on. The most 
interesting question is a possibilty of the TI development 
behind shock waves. The an alysis of stability in the post- 
shock flow was performed b y IShchekinovl (|l978l . Il979l ) and 
IChevalier fc Imamural (|l982l ). Further progress is reached by 
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lYamadafc Nishil |200ll ). who have analyzed the thermal in- 
stability behind radiative shock waves and studied its possi- 
ble role in the fragmentation. In the numerical simulations 
the role of thermal instability behind radiative shock waves 
was studied by many authors (e.g., Binette et al. 1 1 19851 ; 



iDavid et al~l 1 19881 : ISutherland et al. I |2003h . In the above 
metioned papers the radiative shock waves with velocities 
~150 - 1500 km s _1 are considered. On the other hand the 
thermal/thermo-reactive instabilities are analyzed in a dif- 
fuse in terstellar photoionized gas with temperatur e below 
10 4 K jCorbelli fc Ferrara|[l995l ; ISmith et al. 1120081 ). Shock 
waves with velocities ~ 50 — 150 km s _1 are quite impor- 
tant for the evolution of dwarf galaxies. Such shock veloci- 
ties can be assotiated with global star formation processes 
and winds in dwarf galaxies. Thus, it is interesting to study 
both a possibility of the thermal instability and the ioniza- 
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tion and thermal evolution in a gas behind shock waves with 
v ~ 50 - 150 km s _1 . 

Such shock waves are belived to play a significant role 
in the metal enrichment of the intergalactic medium (IGM). 
Metals (heavy elements) produced by stars are transported 
into the IG M by shock waves from galaxies and clusters 
of galaxies JCnedinl ll99Sl; iMadau et al. Il200ll ; IShchekinovl 
120021 ; lAguirre fc Schavdl2009l ; lMeiksinll2009r i. The efficiency 
of metal ejection depends on many factors and parameters, 
e.g. mass of a parent (for metals) galaxy, star formation rate, 
density profile etc. On one hand massive galaxies should pro- 
duce and may throw out a sufficient part of metals because 
of high velocity shock waves. But on the other the escape 
velocity is quite high for such galaxies and metals may be 
confined inside them. Whereas during star formation burst 
dwarf galaxies c an expel major par t of their metal prod- 
ucts in the IGM jFerrara et al. 1120001 s ) . The velocity of shock 
waves produced by dwarfs is about several dozens km s . 
Low-velocity shoc k waves do not lead to the radiative pre- 
cursor formation (|Dopita fc Suther land 199fj). The ionizing 
flux produced by the stellar population of a parent galaxy 
is expected to be insignificant due to the burst character of 
star formation process (the massive part of the stellar popu- 
lation explodes as supernovae in a short timescale, whereas 
low mass stars do not produce sufficient number of ionizing 
photons). Therefore, the ionization and thermal evolution of 
gas behind low-velocity shocks is mainly governed by colli- 
sions between atoms, ions and electrons. 

The dynamics of interstellar/intergalactic gas in gen- 
eral may be understood by studying the ionization states 
of metals. Indeed, the CIV, NV, and OVI ions are sensi- 
tive t racers of hot gas with T about several times 10 s K 
(e.g.. lEdear fc Cheyalierl Il98rf) . Using theoretical models 
llndebetouw fc Shulll |2004al ) have investigated the ioniza- 



tion ratios of the Li-like absorbers C I V, NV and OVI in 
the Galactic halo, I Gnat fc Sternberg! (|2009h have consid- 
ered the ionic column densities as tracers of the thermal 
and ionization evolution behind strong radiative shocks 
with velocities more than hundred km s - . Shock waves 
with diffe rent velocities can be i dentified using the spec- 
tral lines (|Cox fc Raymond! I1985T ). The growth of thermal 
instabilities can be dete cted by the associa ted X-ray and 
optical-ultraviolet lines (|David et al. I Il988h . The analysis 
of the observed column densities of ions can give infor- 
mation ab o ut int e rstellar and intergala c tic structures (e.g., 
Fox et al. I 120051; Islmcoe et al. I 120061; I Gnat fc Sternberg 



20071 ; lAgafonova et al. 1 120071 ; IVasiliev et al. I l2010h . How- 
ever, significant difficulties are connected to adequate in- 
terpretation of the obsevational data. Theoretical models of 
the ionization mechani sms give diffe rent predictions for ions 
and their ratios (e.g.. [Spitzerl Il99r3 ) . and usually the ionic 
ratios do not remove the ambiguities regarding t he ioniza- 
tion conditions (e.g.. llndebetouw fc Shull l2004rj ). But we 
can expect that ionizaton states of metals may help to rec- 
ognize shock waves with v ~ 50 — 150 km s _1 and determine 
physical conditions in postshock gas. 

In this paper we study thermal instability in a colli- 
sionally controlled gas behind shock waves with v ~ 50 — 
150 km s~ . A possible influence from external ionizing ra- 
diation field will be considered elsewhere. The paper is orga- 
nized as follows. In Section 2 we briefly describe the details 



of the model. In Sections 3 and 4 we present our results. In 
Section 5 we summarize our results. 



2 MODEL DESCRIPTION 

Here we briefly describe our method of calculation. The 
full description of the meth od and t he re ferences to the 
atomic data can be found in IVasilievI (|201ll ). We study the 
ionization and thermal evolution of a lagrangian element of 
cooling gas. In our calculations we consider all ionization 
states of the elements H, He, C, N, O, Ne, Mg, Si and Fe. 
We take into account the following major processes in a col- 
lisinal gas: collisional ionization, radiative and dielectronic 
recombination as well as charge transfer in collisions with 
hydrogen and helium atoms and ions. 

The system of time-dependent ionization state equa- 
tions should be complemented by the temperature equation. 
Neglecting the change of number of particles in the system 
(for fully ionized hydrogen and helium it remains approxi- 
mately constant) the gas temperature is determined by 



dT 
~dt 



Anks 



(1) 



where n, n e and uh are the electron and total hydrogen 
number densities, A(xi, T, Z) is cooling rate, A is a constant 
equal to 3/2 for isochoric and 5/2 for isobaric cooling, ks is 
the Boltzman constant. Cooling is isobaric when the cooling 
time is much greater then the dynamical time of a system, 
tc/td 3> 1 and it is isochoric when the time ratio is opposite, 

tc/t d < 1. 

The total cooling rate is calcu lated using the phot oion- 
ization code CLOUDY (ver. 08.00. iFerland et al. Ill998f ). For 
the solar metal li city w e adopt the abundances reported by 
lAsplund et al. I |2005l ), except Ne for whic h the enhanced 
abundance is adopted (|Drake fc Te sta 2005). In all our cal- 
culations we assume the helium mass fraction Yhb = 0.24. 
We solve a set of 96 coupled equations (95 for ionization 
states and one for temperature) using a Variable-coefficien t 
Ordinary Differential Equation solver (|Brown et al. |[l989l ). 

We consider pure collisional ionization model, so we 
should constrain shock wave velocity by a value that does 
not lead to the radiative precursor formation. The precusor 
is photoionized by the radiation field emitted by a shocked 
gas, and consequently the initial ionization states in gas be- 
hind a shock front depends on the spectrum emitted by 
a shock wave. A stable photoionization precursor will be 
formed when the ionization front velocity in the gas ap- 
proaching the shock becomes larger than the shock veloc- 
ity. This conditio n is satisfied for shock veloc ity higher than 
«* £ 175 km s" 1 (|Dopita fc Sutherland ^d - ) or shock tem- 
perature T, = 3m p v*/16k > 7 x 10 5 K. Thus, we study the 
shock waves with temperature T 3 ^ T». We constrain our 
calculations by lower value v — 155 km s _1 corresponded to 
T s ~ 6 x 10 5 K. We assume that the gas temperature just be- 
hind a shock front instantaneously becomes equal to T 3 , but 
the ionization composition relaxes from that corresponded 
to temperature T; before the front. During the relaxation 
period a gas cools more efficiently due to higher collisional 
ionization rates of low ionic states existed in the postshock 
gas, and the ionization composition tends to fit such ther- 
mal evolution. In general, for fixed other paramters (den- 
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Figure 1. The hydrogen and helium ionization states for a 
gas isochorically cooled with 10 -3 , 0.1 and solar metallicities are 
shown by dot, solid and dash lines, respectively. 




log T, K 



Figure 2. The contributions to the total cooling rate from each 
chemical element for T; = 8 X 10 4 K (upper panel) and T; = 
3 X 10 5 K (lower panel) in a gas with solar metallicity. The cooling 
rate for the fiducial model, T; = 10 s K, is shown by thin dotted 
line. 



COOLING RATES AND THERMAL 
INSTABILITY 



sity, metallicity) the ionization composition tends to that at 
a corresponding temeperature in a gas collisionally cooled 
from T = 10 8 K. However, depending on the relation be- 
tween Ti and T s the relaxation period needs different time. 

The ionization times of H and He are longer than these 
of metals, so that the influence of H and He on the thermal 
evolution of gas just behind a shock is expected to be domi- 
nant. As an example Figure [T] presents the H and He ioniza- 
tion states in a gas isochorically cooled from T = 10 8 K for 
1CP 3 , 0.1 and solar metallicities. There is no significant de- 
pendence of the hydrogen states on metallicty, in contrast to 
that the helium states show remarkable distinction. The con- 
siderable dependence of the H and He states on temperature 
is clearly seen. Thus, taking different ionization composition 
corresponded T = T in a gas before shock wave one can ex- 
pect differences in the thermal evolution of the postshock 
gas. But the closer T is to T s , the smaller distinction of the 
thermal evolution is from that for a gas collisionally cooled 
from T = 10 s K. 

A gas temperature before shock wave can vary in wide 
range. For example, a gas can be relaxed from some previous 
shock interaction or heated by galactic or extragalactic ioniz- 
ing radiation. So it is difficult to choose more or less common 
initial temperature. The gaseous temperature in HII regions 
is ~ 2 — 4 x 10 4 K, the typical temperature of the IGM 
is ~ 2 x 10 4 K. To coinstrain calculations we suppose that 
the initial temperature of a gas before shock wave equals 
T = 2 x 10 4 K, and a gas has the ionization composition 
obtained in the non-equilibrium (time-dependent) model at 
this temperature and given metallicity. But in several sets 
of models we study the depedence of the gas evolution on 
T < T s . For these models the ionization composition is cor- 
responded to T. 



The thermal instabili ty criteria for a ga s cooled isochorically 
and isobarically are (|Shchekinovll 197Sl ) 



(2) 



d In T d In T 

respectively. The isochoric and isobaric cooling rates for 
a gas cooled from T — 10 s K are close to each other 
(|Gnat fc Sternberg! 120071 ). So that a gas is unstable to 
isochoric/isobaric perturbations in very close temperature 
ranges. However, the cooling rate of a gas behind a shock 
may depend on physical condition in a gas before shock 
front, e.g. ionization composition and temperature of a gas, 
so that the temperature range, where a gas becomes ther- 
mally unstable, can be modified. 

Figure [3] presents the isochoric cooling rates, thermal 
instability criterion, /, and characteristic scale for a gas 
with 10 -3 , 0.1 and 1 Zq behind a shock with T s . Below 
we present our calculations in terms of a shock temperature 
T s = 3m p Vs/16k ~ 2.3 x 10 3 vf , where v w = v a /W km s -1 . 
Also for a given metallicity the isochoric cooling rate for 
T — T s — 10 s K, a fiducial m odel, which corresponds to 
the tabulated cooling ra tes (e.g.. lSutherland fc Dopitalll993l ; 
iGnat fc Sternb erg 2007), is plotted by thin solid line. First 
of all we should note that no significant difference between 
isobaric and isochoric cooling rates for T B 6 x 10 5 K is 
found. So here we present figures for the isochoric case only, 
but due to the different thermal instability criterion we an- 
alyze both cases. 

Figure [3{ upper panels) shows that the cooling rates 
for the whole shock temperature range considered here, 
T s = 0.5 — 6 x 10 5 K, tend to the fiducial one and equal 
to it at low temperature. For gas with Z = 10 -3 Zq this oc- 
curs almost at the T s due to fast ionization of hydrogen for 
T s ^ 7 x 10 4 K and both H and He ionization for higher ini- 




Figure 3. The cooling rates (upper panels), thermal instability criterion (middle) and characteristic scale (lower) for gas wih 10 3 , 0.1 
and 1 Zq metallicities (from left to right). Upper panels. The cooling rates for T s = 5 X 10 4 , 6 X 10 4 , 7 X 10 4 , 8 X 10 4 , 9 X 10 4 , 10 5 , 
2 x 10 s , 3 x 10 5 K are depicted by solid lines from left to right, respectively. For solar metallicity (right column of panels) the cooling 
rates for T a = 4 x 10 5 , 5 x 10 6 and 6 X 10 5 K are added. The cooling rate for the fiducial model is shown by dash line. Middle panels. 
The same for the thermal instability criterion. Two horizontal lines corrcspod to the criterion for isobaric (upper) and isochoric (lower 
line) perturbations. A gas is unstable below the lines. Lower panels. The column density of a thermally unstable gas layer, Nt = n/X, 
for isobaric (shown by thin solid lines) and isochoric (shown by thick solid lines) perturbations, where At = c s t coo i is the thermal length. 
The lines from bottom to top correspond to higher shock temperature, T s . For Z = 10 — 3 Zq the lines coincide. 



tial temperature (the vertical part of lines in panel a). For 
higher metallicity other chemical elements, mainly carbon 
and oxygen, can dominate in cooling. This is clearly seen 
in Figure [2l which presents the contributions to the total 
cooling rate from each chemical element for T s = 8 X 10 4 K 
(upper panel) and T s = 3 x 10 s K (lower panel) in a gas with 
solar metallicity. One note that the difference between cool- 
ing rates behind shock waves and the rate for the fiducial 
model reaches maximum for solar metallicity and depends 
strongly on T s , but for any T 3 the cooling rates coincide at 
T < 3 x TO 4 K. Thus, we can conclude that the cooling rate 
in a gas behind a shock wave differs considerably from the 
cooling rate for the fiducial model, T s = 10 8 K, which is 
usually used to study the thermal evolution of a gas. 

Panels (b) present the / value for the considered mod- 
els. Two horizontal lines correspond to the isobaric/isochoric 



thermal instability criteria. Note that both criteria are sat- 
isfied for the fiducial model at T £ 3 x 10 4 K for any metal- 
licity considered here. For Z = 10~ 3 Zq both isobaric and 
isochoric criteria for the models with T s — 0.5 — 6 x 10° K are 
also satisfied (the region below horizontal lines) in the same 
temperature range as that for the fiducial model, except for 
very early evolution of a gas behind shock wave (almost ver- 
tical tails of the cooling curves correpond to that period of 
the evolution). Increase of metallicity leads to widening tem- 
perature range, where the isochoric criterion is not satisfied, 
but does not change the range for the isobaric case. A gas 
evolved from T s = 0.5 — 6 x 10 5 K is thermally unstable to 
isobaric perturbations during almost full its evolution. In a 
gas cooled isochorically the value / for the solar metallicity 
demonstrates a significant difference from that for the fidu- 
cial model. The temperature range, where a gas is unstable 



Thermal instability in the collisionally cooled gas 5 



1.5 



1 - 




-1.5 1 1 1 1 1 1 1 1 ' — 

4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 

log T, K 



Figure 4. The critical column density, Dt r n, for the models with 
T s = 5 x 10 4 , 6 x 10 4 , 7 x 10 4 , 8 x 10 4 , 9 x 10 4 , 10 5 , 2 x 10 5 , 
3 X 10 5 , 4 X 10 5 , 5 X 10 6 and 6 X 10 5 K (thick dashed lines from 
left to right, respectively) and for the fiducicial model (thin solid 
line) at solar metallicity. 



behind shock waves with T s > 6 x 10 s K, is very close to 
that for the fiducial model. But such temperature range be- 
comes narrower for lower T s , e.g. for T s = 3 x 10 5 K a gas 
is unstable at T ~ 2.6 - 5.4 x 10 4 K, and for T s = 10 5 K 
it proceeds at T ~ 2.8 - 4 x 10 4 K. For T s < 10 5 K a gas 
remains thermally stable to isochoric perturbations during 
full its evolution. 

Panels (c) present the column density of a gas layer, 
which can be thermally unstable, Nt = n/Xt, where A t = 
Cstcooi is the thermal length. We have plotted the depen- 
dence of Nt, where the instability criterion is satisfied, 
namely, for any T s at T £ 2 x 10 4 K in the isobaric case 
(shown by thin solid lines) as well as for T s £ 10 s K in the 
isochoric one (shown by thick solid lines). Note that the cri- 
terion for the fiducial model with T, = 10 8 K is satisfied for 
r^2x 10 4 K (dash line). The difference between column 
densitites for T s = 0.5 — 3 x 10 s K and these for the fiducial 
model reaches a factor of 2-3 for Z ^ 0.1, it is maximum 
for early stages of the evolution, where the cooling rate is 
changed considerably. 

Following Gnat & Sternberg (2007) we consider the 
transition from isobaric to isochoric cooling for the inter- 
action a shock wave with a cloud in terms of the critical 
column density, D tr n, and temperature, where Dt r is a crit- 
ical size of cloud. Figure [4] present such dependence for the 
set of models with T s = 0.5 — 6 x 10 5 K (solid lines) and 
for the fiducial model (dash line) for solar metallicity. The 
cooling becomes isochoric for temperatures and cloud col- 
umn densities left and above of the curves. Here D tr n oc c s , 
so that Dt r n should be multiplied by the Mach number M 
for a shock wave with v s = Mc s - 

In the above-considered set of models we start from the 
ionization composition corresponded to the Ti = 2 x 10 4 K. 
But a gas behind shock wave can have another ioniza- 
tion composition. Figure [5] presents the isochoric cooling 
rates, thermal instability criterion and characteristic scale 
for T s = 10 K depending on the initial ionization com- 
position, which is corresponded to Tj = (2 — 9) x 10 4 K. 
Figure [5j upper panel) shows that the cooling rates gradu- 
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Figure 5. The same as in Figure \3\ but for the shock wave 
with T s = 10 5 K depending on the initial ionization composition, 
which is corresponded to T t = 2 x 10 4 , 3 X 10 4 , 4 X 10 4 , 5 X 10 4 , 
6xl0 4 , 7xl0 4 , 8xl0 4 , 9xl0 4 K from top to bottom, respectively. 

ally approach to the fiducial one, the difference is negligible 
for Ti ^ 7 x 10 4 K. The thermal instability criterion and the 
value nXt also demonstrate the same behavior. 



4 IONIZATION STATES AND THEIR RATIOS 

Figures [6]|8] show the A'nv/^ovi, Nciv /Nshv and 
Nciv/Noyi ionization ratios in a gas with 10~ 3 , 0.1 and 
IZq metallicities behind a shock wave with T s . The tem- 
perature before a shock wave is taken Ti = 2 x 10 4 K. Here 
we consider isochoric ionization ratios, although the differ- 
ence from isobaric ratios is non-negligible, their temperature 
dependences show similar behavior. 
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Figure 6. The iV NV /./Vovi> Ncrv/Nsnv an d A r civ/A r ovi ioniza- 
tion ratios for 10~ 3 Zq metallicity gas with Ti = 2 X 10 4 K behind 
a shock wave with T s = 5 X 10 4 K (thin solid line), 7 X 10 4 K 
(dash line), 9 X 10 4 K (short dash line), 10 5 K (dot line), 2 X 10 5 K 
(dash-dot line), 3 X 10 5 K (short dash-dot line), 4 X 10 5 K (short 
dash-dash line), 5 X 10 K (dash-dot-dot line) and the fiducial 
model, T a = 10 s K, (thick solid line). 



The ionization ratios for the metallicities considered 
here demonstrate a significant dependence on T s and dif- 
ference from the ratios in the fiducial model. The differ- 
ence in some temperature ranges reaches several orders of 
magnitude. The difference of the Nciv/Novi ratio increases 
with the growth of T„ and reach maximum (up to two or- 
ders of mag nitude) for T 3 = 2 - 3 x 10 5 K in a gas with 
Z — 1CP' ! — IZq (see the lower panels in Figures [6][8j) . Fur- 
ther increase of T s diminishes the difference rapidly, and the 
ratios for T s = 4 — 5 x 10 J K almost coincide with the fiducial 
one. Stronger shock wave (higher T 3 ) leads to the ionization 
of higher states, and the difference from the fiducial model 
should decrease with the growth T„. Two other ratios exhibit 
similar behavior. 

The difference of the ionic ratios from the fiducial ones 
is maximum for the lowest metallicity, Z = 10~ 3 Zq. For ex- 
ample, the iVcrv /Novi ratio for T s = 2 x 10 5 K differs about 
four orders of magnitude from that in the fiducial model. In 
a solar metallicity gas this difference reaches two orders only. 
Obviously, the ionic fractions decrea se fastly in low metallic- 
ity gas (e.g. lGnat fc Sternberg|2007l ). and the recombination 
lag increases with the metallicity growth. So the compari- 
son of two ions with large gap between ionization potentials, 
like CIV and OVI (47.9 and 113.9 eV for CIII^CIV and 
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Figure 7. The same as in Figure[6] but for Z = O.IZq. 



OV— s-OVI, respectively), gives higher values for lower metal- 
licity. Whereas the ratio between CIV and SilV demonstrate 
weaker dependence on metallicity due to proximity of their 
ionization potentials (34 eV for Siffl— >SiIV) . 

Such a strong dependence on T s may lead to some un- 
certainly or inaccuracy of determining the physical con- 
ditions in the postshock gas using the tabulated colli- 
sional ionization states (e.g. ISutherland fc Dopital 1 19931 : 
I Gnat fc Sternberdl2007h . In the solar metallicity gas the 
most remarkable dependence is within the temperature 
range, where the postshock gas can be thermally stable for 
isochoric perturbations (see Figure [3j and the transition 
from isobaric to isochoric cooling is more probable (Fig- 
ure Moreover, the postshock material can be efficiently 
mixed due to the hydrodynamical instab ilities (|Slavin et al. I 
1 19931 : Ide Avillez fc Breitschwerdt 120091 ) , and a degree of un- 
certainly of determining the physical conditions may in- 
crease. So that the strong dependence of ionic states can 
hardly help to dicriminate this transition, and it should be 
taken into account for correct interpretation of observational 
data and synthetic data models obtained from numerical 
simulations. 

Also we have analyzed the dependence on the initial 
ionization composition of a gas. Figure [5] shows the same 
ionization ratios, AW/A'ovi, Nciv/Nsuv and A r civ/A^ovi, 
in a gas behind a shock wave with T s = 3 x 10 5 K for differ- 
ent initial ionization composition corresponded to Ti. The 
increase of Ti leads to less difference of the ionic ratios from 
these in the fiducial model. As it is mentioned above such 
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Figure 8. The same as in Figure [6] but for Z = Zq. 



difference is more pronounced for ions with larger deviation 
in their ionization potentials. 

Figure [10] (lower panel) presents the dependence of 
Nciv / Novi versus Nnv/Novi column densities on the tem- 
perature behind a shock wave with T 3 for a gas with so- 
lar metallicity and the initial ionization composition corre- 
sponded to Ti = 2 x 10 4 K. All tracks start from the same 
point corresponded to the initial ratio at T = 2 x 10 4 K. 
Further evolution depends on collisional ionization rate of 
ions behind a shock with T s . As it is expected from the 
above-mentioned results the tracks differ significantly from 
the fiducial model: the largest deviation can be found for 
T 3 ~ 2 - 3 x 10 5 K, while the tracks for T s < 10 5 K are close 
to each other, and the tracks for T 3 ^ 4 x 10 K becomes 
closer to that for the fiducial model. The difference becomes 
negligible at T s > 6 x 10 s K. 

Finally we again consider the dependence on the initial 
ionization composition of a gas. Figure \W\ (upper panel) 
shows the dependence of the column densities behind a 
shock wave with T s = 3 x 10 5 K on the initial ionization 
composition corresponded to T. We should note that the 
tracks for the shock models tends gradually to the fiducial 
track, but these tracks coincide at T 4 x 10 4 K. How- 
ever, at T ^ 4 x 10 4 K a significant difference can be found 
even for the track T = 10 5 K. So that even small deviation 
of the initial ionic composition leads to remarkably differ- 
ent tracks. In both panels we add the observational data 
for t he high velocity clouds in the Galactic halo (see Table 
8 m iFox et al. I l2005h . Although several tracks are within 
the errorbars of the observational data, the majority of the 
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Figure 9. The N N y/Novi, Nciv/Nsnv and Nciv/Novi ion- 
ization ratios for solar metallicity gas with Ti = 2 X 10 4 K (thin 
solid line), 3 X 10 4 K (dash line), 5 X 10 4 K (short dash line), 
7 X 10 4 K (dot line), 9 X 10 4 K (dash-dot line), 10 5 K (short 
dash-dot line) behind a shock wave with T s = 3 X 10 5 K. The 
ratios for the fiducial model, T a = 10 8 K, are depicted by thick 
solid line). 



points cannot be fitted by any track presented in Figure [TOl 
These observational points are belived to be explained by 
turbulent mixing behind shock w aves (SI avin et al. I Il993l ; 
Ide Avillez fc Breitschwerdt Il2009h . 



5 CONCLUSIONS 

In this paper we have studied the thermal instabil- 
ity, nonequilibrium cooling rates and ionization states in a 
collisionally controlled gas behind shock waves with v ~ 
50 - 150 km s" 1 (T s ~0.5-6x 10 5 K). Such shock waves 
do not lead to the radiative precursor formation, and gas 
evolution is governed by collisions only. 

Our results can be summarized as follows. 

• The cooling rate in a gas behind shock waves with T s ~ 
0.5 — 6 x 10 5 K differs considerably from the cooling rate for 
a gas cooled from T — 10 8 K. 

• The temperature range, where the postshock gas is 
thermally unstable, is significantly modified and depends on 
gas metallicity. For Z ~ 10~ 3 Zq both isobaric and iso- 
choric criteria for shock waves with T s = 0.5 - 6 x 10 5 K 
and T s = 10 8 K are satisfied in the same temperature 
range. Increase of metallicity leads to widening tempera- 
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logf(NV/OVI) 

Figure 10. Lower panel. The dependence Nqiv/Novi versus 
Nnv/Noyi for solar metallicity gas with Ti = 2 X 10 4 K behind 
a shock wave with T s = 5 X 10 4 K (label '54'), 9 X 10 4 K (label 
'94'), 10 5 K (label '15'), 2x 10 5 K (label 25), 3x 10 5 K (label '35') 
and the fiducial model, T s = 10 s K, (label 'cc'). Upper panel. The 
dependence Nciv /Novi versus TVnv /Novi for solar metallicity 
gas with Ti = 2xl0 4 K (label '24'), 3xl0 4 K (label '34'), 5xl0 4 K 
(label '54'), 7 x 10 4 K (label '74'), 9 x 10 4 K (label '94'), 10 s K 
(label '15') behind a shock wave with T s = 3 x 10 5 K. The ratios 
for the fiducial model, T s = 10 s K, are marked by label 'cc'. Gas 
temperature is indicated by gray scale along the trajectories: from 
hot (light gray) to cold (black) gas. The data poi nts show the ioni c 
ratios observed in metal absorbers (see Table 8 iFox et al. 112005]) . 

ture range, where the criterion is not fulfilled. For the solar 
metallicity a gas behind a shock with v ~ 50 — 120 km s _1 
(T s = 0.5 — 3 x 10 5 K) is thermally unstable to isobaric per- 
turbations during almost full its evolution. But the temper- 
ature range, where a gas is unstable to isochoric perturba- 
tions, becomes narrower for lower T B , and a gas remains ther- 
mally stable to isochoric perturbations behind shock waves 
with v 65 km s _1 (T B 10 5 K) during full its evolution. 

• The column density of a gas layer, which can be ther- 
mally unstable, also depends on gas metallicity. The differ- 
ence between column densitites for T s = 0.5 — 3 x 10 5 K and 
these for a gas cooled from T = 10 8 K reaches a factor of 
2-3 for Z > 0.1Z Q . 

• The transition from isobaric to isochoric cooling for a 
gas with Z ^ 0.1 Zq behind shock waves with T s = 0.5 — 
3 x 10 5 K proceeds earlier (at lower column density layer) 
than that occurs in a gas cooled from T = 10 s K. 

• The ionic ratios in a gas with Z ~ 10~ 3 — IZ@ behind 
shock waves with T s £ 4 x 10 s K demonstrate a significant 
dependence on T s and difference from the ratios in a gas 
cooled from T = 10 8 K. The difference becomes negligible 
at T s ^ 6 x 10 5 K in the metallicity range considered here. 
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